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ABSTRACT 

We calculate the non-linear matter power spectrum using the 3rd-order perturbation theory without 
ignoring the pressure gradient term. We consider a semi-realistic system consisting of two matter 
components with and without pressure, and both are expanded into the 3rd order in perturbations 
in a self-consistent manner, for the first time. While the pressured component may be identified with 
baryons or neutrinos, in this paper we mainly explore the physics of the non- linear pressure effect 
using a toy model in which the Jeans length does not depend on time, i.e., the sound speed decreases 
as a""'^/^, where a is the scale factor. The linear analysis shows that the power spectrum below the 
so-called filtering scale is suppressed relative to the power spectrum of the cold dark matter. Our 
non-linear calculation shows that the actual filtering scale for a given sound speed is smaller than the 
linear filtering scale by a factor depending on the redshift and the Jeans length. A ~ 40% change is 
common, and our results suggest that, when applied to baryons, the temperature of the Inter-galactic 
Medium inferred from the filtering scale observed in the flux power spectrum of Lyman-a forests 
would be underestimated by a factor of two, if one used the linear filtering scale to interpret the data. 
The filtering mass, which is proportional to the filtering scale cubed, can also be significantly smaller 
than the linear theory prediction especially at low redshift, where the actual filtering mass can be 
smaller than the linear prediction by a factor of three. Finally, when applied to neutrinos, we find that 
neutrino perturbations deviate significantly from linear perturbations even below the free-streaming 
scales, and thus neutrinos cannot be treated as linear perturbations. 
Subject headings: cosmology : theory — large-scale structure of universe 



1. INTRODUCTION 

Pressure plays an important role for the structure for- 
mation in the universe. Pressure determines the Jeans 
scale, A J, below which the growth of structure slows 
down, and eventually stops and oscillates: while fiuctua- 
tions in the cold dark matter (CDM) and the pressured 
component evolve in the same way above the Jeans scale, 
their evolutions are significantly different below the Jeans 
scale. 

The dominant source of gravity is CDM, which is cold 
and its velocity dispersion is negligible before the collapse 
of halos. However, the sub-dominant matter components 
- baryons and neutrinos - have significant velocity disper- 
sions, which should be included in the calculation when 
precision is required. While the accurate calculations 
have been done for the linear perturbations, the effects 
of the pressure on the non-linear evolution of matter fluc- 
tuations on cosmological scales (~ 10 — 100 Mpc) have 
not been studied very much in the literature. 

We address this issue by calculating the non-linear 
matter power spec trum using the 3rd -order perturbation 
theory (3PT; see iBernardeau et al.ll2002l for a review), 
with the pressure gradient term in the Euler equation ex- 
plicitly included. This enables us to study the effects of 
the pressure on the non-linear evolution of matter fluc- 
tuations in a self-consistent manner. 

The rest of this paper is organized as follows. In §[2l we 
find the linear, second-order, and third-order solutions 
of the coupled continuity, Euler, and Poisson equations 
for two matter components with and without the pres- 
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sure gradient. In §[3l we calculate the non-linear matter 
power spectrum from the solutions obtained in § [21 In 
§ m we compare ou r full 3PT calcula tion with the ap- 
proximation used by lSaito et al.l ()2008[ ) for the effects of 
massive neutrinos on the matter power spectrum. Fi- 
nally, in § we discuss the implications of our results 
for a few practical astrophysical and cosmological appli- 
cations. In Appendices we give the detailed derivations 
of the 3PT results used in the main body of the paper. 

2. THIRD-ORDER PERTURBATION THEORY WITH 
PRESSURE 

2.1. Basic Equations 

The main goal of this paper is to find the perturbative 
solutions for the CDM density contrast, 5c, for which the 
pressure gradient is ignored, and the density contrast of 
another matter component, 6b, for which the pressure 
gradient is retained. This component may be identified 
with baryons (hence the subscript "b" ) or neutrinos, de- 
pending on the sound speed one uses in the Euler equa- 
tion.^ 

The equations that we are going to solve include two 
continuity equations: 

,5,(x,r)+V - [(1 + ^c(x,t))v,(x,t)] = 0, (1) 
,5b(x, r) + V • [(1 + <5fc(x, r))vb(x, r)] = 0, (2) 

^ While we use "b" to denote the pressured matter component 
throughout this paper, we do not always mean baryons, but we we 
always refer to a general matter component with pressure. 
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two Euler equations: 
Vc(x, r) + [Vc(x, r) • V]vc(x, r) 



Vb(x, t) + [vb(x, t) • V]v6(x, t) = Vb{x, t) - V0(x, r) 



Cs(x, r)VJb(x,T) 
1 + (56(x,t) 



,(4) 



and one Poisson equation: 

V20(x,r) =4^Ga2[p,(T)J,(x,T) +p;,(r)5;,(x,T)],(5) 

where 6i = (pi — pi)/ pi is the density contrast of a matter 
component i — {c,b), p the background matter density, 
a the scale factor, the pecuhar velocity field of a mat- 
ter component i, (f> the gravitational potential, and Cs 
the sound speed of the matter component with pressure. 
Here, the dots denote the partial derivatives with respect 
to the conformal time, r, i.e., 5 = 85 /dr, and V denotes 
the partial derivatives with respect to the comoving co- 
ordinates. 
We rewrite the Poisson equation as 



6 



(6) 



where we have assumed an Einstein-de Sitter (EdS) uni- 
verse (we shall generalize the results to other cosmolog- 
ical models later), for which the energy density of the 
universe is dominated entirely by the matter density, and 
a oc . The background Fricdmann equation is given by 



(7) 



We have also defined the total matter fluctuation, 6, 
which is given by 



J(x,t) 



^ Pc(T)^e(x, r) + pb{T)Sb(pC, t) 
~ Pc{T)+pb{T) 



where fc = Pc/{Pc + Pb) = ^d^m, and fb = Pb/{Pc + 
pb) — flb/^m- For an EdS universe, flm — 1. 

Taking the divergence of the Euler equations, we ob- 
tain the equations for the velocity divergence fields, 
9i = V • Vi. Moving non- linear terms to the right hand 
side (RHS) of the equations and using the Poisson equa- 
tion, we obtain 



jc(x, r) + edx, t) = -V • [5c(x, t)vc(x, t)], 

(56(x, r) + 9bix, t) = -V • [(56(x, r)vb(x, t)], 
2 6 

T 

= -V-{[v,(x,r) • V]v,(x, r)}. 



= -V-{[vfc(x,T) • V]vb(x,T)}-V- 



l + (5h(x, r) 



(9) 
(10) 

(11) 



Note that the second term in the RHS of eq. [12] still 
contains the linear order term. All the other terms in 
the RHS of the above equations are non-linear. 

We shall simplify the pressure term, the second term 
in the RHS of eq. [12], as follows. First, we shall assume 
that the sound speed is homogeneous, i.e., Vc^ — 0. See 
iNaoz fc Barkanal (|2005f) for the analysis of linear pertur- 
bations with VCg ^ 0. Second, we expand the pressure 
term to the 3rd order in perturbations: 



Pb 



ySb 

1 + 6 



V5b - 5bV5b + SlVSb + OiSt).{13) 



/c<5c(x,r)+/b5b(x, 
(8) 



r), 



Going to Fourier space, we obtain 



k{k,T) + eb{k,T)^- 



1 



(27r)3 
1 

(2^ 



dqidq2^D(qi + q2 - k) — —9c{qi,T)Sc{q2,T), 

1i 

dqi dq2 (5n (qi + q2 - k) — — 6*6 (qi , r) (5b (q2 , t) , 

1i 



i{k,T) + -9,{k,T) + -^~S{k,T) 
T 

9b{k,T) + -eb{k,T) + -^d{k,T) 



1 



(27r)3 
1 

(2^ 



dqidq2(5D(qi + q2 - k) 



dqi(iq2(5n(qi -f qj - k) 



fc^(qi ■ q2) 
fe^(qi ■q2) ^ 



^c(qi,T)6'c(q2,T) 



(qi,r)6'b(q2,T) 



where 



cUT)V5b{x,T) 

l + (5b(x, r) 



(k) 



~SbO^) 



C^,iT)V6bix,T) 

l + (5b(x, r) 
1 



(k), 



2(27r) 



(14) 
(15) 
(16) 

(17) 



dqi(iq25b(qi,r)Jb(q2,T)(5r,(qi +q2-k) 



+ 



3(27r)e 



dqidq2dq3(5b(qi,r)(5b(q2,T)^b(q3,T)fo(qi + q2 + qa - k) 
I ' 



(18) 



In the subsequent subsections we shall solve these cou- 



pled equations perturbatively. Hereafter we shall omit 
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the tildes on the perturbation variables in Fourier space. 



2.2. Linear Order Solution: Jeans Filtering Scale 
In the linear order, one finds 



^i,b(k,T) 

U2 2 



?i,c(k,r) =0, 



(?i.b(k,' 
--^i,c(k,r) 

T 

--^i.b(k,T) 



0, 



fc^c^(T)5i,b(k,r) =0, 



(19) 
(20) 

(21) 
(22) 



where the subscripts "1" mean that these quantities de- 



note the first-order perturbations, and 5i 
fbSi^b- We rewrite eq. [22] as 



fcSl 



9i,6(k, T) + -ei,bCk, t) 

T 



6 



-'^i,b(k, 



9i,6(k, r) -I- -0i.b(k, r) 

T 



= 0, 



(23) 



where we have used the usual definition of the Jeans 
wavenumber, kj: 



kj{r) ^ 



V6 



(24) 



Cs{t)t' 

The Jeans wavenumber divides the solutions for Si j, into 
two classes: the growing solution for k ^ kj, and the 
oscillatory solution for k ^ kj, when there is no CDM, 
i.e., fb = ^ and i5i = Si^b- When i5i ^ 6i^b, the Jeans 
wavenumber does not provide a dividing scale for the 
solutions of Si^b- 
The Jeans wavenumber depends on the temperature of 



the matter component "b" as kj oc T, 



-1/2, 



'T thus, kj 
depends on time in general, kj — kj{T). However, in or- 
der to simplify the problem and obtain physical insights 
into the effects of pressure on the non-linear growth of 
structure, we shall assume that kj is independent of time, 
which requires that the matter temperature evolve as if 
the matter were coupled to radiation, T;, oc 1/a oc l/r^. 
This is not a realistic assumption especially in a low red- 
shift universe where baryons are decoupled from the ra- 
diation background and neutrinos are non-relativistic - in 
both cases the temperature evolves as Tf, oc 1 /a^ oc l/r"* 
and thus kj evolves as fcj oc r oc a^/^, for the adiabatic 
evolution. 

We shall solve the above coupled linear equations iter- 
atively: as CDM is always the most dominant source of 
gravity, the zeroth-order iterative solution may be found 
by setting Si Si^c (i-e., /c 1). We find the solution 
for the ratio of the density cont rasts, which is often c alled 
the "Jeans filtering function," (jGnedin fc HuilfT998f l 



^ _ '5l,b(k,T) 

5i(k,T) = 



<5i,c(k,T) 



(25) 



which should be a decreasing function of k due to the 
effect of pressure. At the zeroth-order of iteration, the 



CDM density contrast grows as 

df'l (k, t) oc a oc t^, 
and thus the equation for gi simplifies to 



(26) 



'^r(k,r) + ^(l + |,),r(k,.) 



(27) 



The solution for gi (k, r) must be normalized such that 
gi{k,T) — > 1 as fc ^ 0. We find 



9f\k,r) = 



1 



O 



where 



24 / fc2 



(28) 



(29) 



The second term is a decaying mode, whose amplitude is 
set by the initial condition, e.g., at the epoch when the 
baryon temperature was raised (by, say, cosmic reion- 
ization) to the point where the pressure became impor- 
tant, or at the epoch when the neutrinos became non- 
relativistic. 

Ignoring the decaying mode (although we shall come 
back to this later), we have the zeroth-order solution: 



9f\k)^ 



1 



fc2 • 

IT 



(30) 



At the first-order iteration we have the pressure feed- 
back on the growth of CDM. The evolution of s'^], de- 
pends on fc, and is given by 



5i;^(fc,r)(xr"W, 



(31) 



where 



n(fc) = - 



-l±5^1-|/,(l-5r^(fc)) 

2-|A[l-5f^(fc)] 
-3+|/,[l-5j°Hfc)] " 



(32) 



The second equality is valid for /^[l — ^^''^(fc)] <C 1. The 
growing mode solution is given by 



n+{k) 



6 



fb[l~9f\k)]. 



(33) 



As g''^\k) 1 and for fc — > and oo, respectively, the 
large-scale and small-scale limits of the growing mode 
solution is (see, e.g.. Sec. 8.3 of IWeinberi l2008l . for a 
recent review) 



^^il+ik, r) oc r'^ oc a, 
5^il+{k,T)oiT'^^i^>' oca^"^-^" 



fc < fcj, (34) 
fc>fcj. (35) 



The growth of 5i c on the spatial scales below the Jeans 
scale is suppressed relative to that of the large-scale 
modes. 
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Taking the first order iteration solution for s[^l^ into 
account, tlie first order iteration equation for g[^^ is 



9['\k,r) + - 



(1)/, N 6(1-/6) 

5i (k,r) = 2 , 

T 

(36) 

whose growing mode solution (with the normalization 
that g[^^ ^ 1 for fc ^ 0) is 



fe2 

IT 



M2^gf'\k)] 

i-h 



1 - + § 



1 - 



(37) 



This iteration converges quickly for < 0.5, and fur- 
ther iterations are not necessary. The largest difference 

between gf'\k) and g[^\k) occurs as k/kj oo, and 
is 100% for /f, = 0.5. If the component "b" is identified 
with baryons, /b ~ 1/6, and the difference is reduced to 

20%. The difference between g[^\k) and g^\k) oc- 
curs at fc ~ kj, and is ~ 4% for /b = 0.5, and 0.2% for 
/{, ~ 1/6. The difference is much smaller for neutrinos. 

To simplify the subsequent analysis, we shall adopt the 
zeroth-order iterative solution for the filtering function, 

gf'^ = 1/(1-1- k'^/kj), and the first-order iterative solu- 
tion for the CDM growth factor, eq. as the solution 
at the first-order in perturbations. This solution is suffi- 
ciently accurate for our obtaining the physical insights. 

Let us comment on the decaying mode that we have 
ignored in obtaining eq. [30j . This decaying mode is an 
oscillatory function at k/kj > 1/(2^/6) ~ 0.2, represent- 
ing the acous tic oscillation of the pressured component 

usse 3 [2000l) . While this ter m is a decaying mode, it 
decays slowly, and is not quite negligible even at low red- 
shift. We show the decaying mode at the zeroth-order 
iterative solution in Fig. [1] 



Ag^^\k,T)^g^°\k,T)- 



1 



fc2 



(38) 



assuming that the pressure became important at = 
10. This figure shows that the decaying mode remains 
important even until z 0; thus, technically speaking, 
ignoring the decaying mode results in an inaccurate form 
of the filtering function. Nevertheless, we shall ignore it 
and adopt gi{k) = 1/(1 + k'^/kj). 

The exact form of gi{k,T) is not so important for our 
purposes. The main goal of this paper is to study how 
non-linearities affect this function. In other words, we 
are interested in how the higher-order filtering functions, 
g„(k, t), are related to the linear one, gi{k, r). One may 
use any forms of gi{k, r) for a better accuracy, depending 
on the problem (baryons or neutrinos). 

2.3. Second and Third Order Solutions 

For the higher order (n-th order) density perturbations 
and velocity-divergence fields, we define the Jeans filter- 




FlG. 1. — Decaying mode solution for the linear filtering function 

at the zeroth-order iteration (/c — > 1), Ag^-^\k,T) = g^-^\k,T) — 

l/(l + fe'^/fcj), where (k, r) is the numerical solution of eq. I I27I I. 

with the initial conditions given by g^^ (k, r, ) = 1 and g^^ (fc, t» ) = 
where t» is the conformal time at z* = 10. The top and bottom 
lines at k/kj ^ 1 are at z = 8 and 0, respectively, and the other 
lines correspond to the intermediate redshifts. 



ing functions such that 

gn(k,T) 
/l„(k,T) 



^ (5n^b(k, r) 
(5„,c(k, r) ' 
^ ^^^^(k, r) 
~ 6'„,c(k,T) ' 



(39) 
(40) 



Assuming that CDM dominates the gravitational po- 
tential, we find the zeroth-order iteration ansatz in an 
EdS universe: 

oo 

5t{k, r) = ^ a"(r)5„,c(k)5„(k, r), (41) 

oo 

Obik, r) = ^ a(T)a"-i(T)0„,,(k)/i„(k, r). (42) 

Tl=l 

Detailed derivations of the non-linear filtering func- 
tions at the second order, (72(k, r), and the third order, 
(73(k, t), are given in Appendix § |B] The second-order 
solution is 



52 (k, r) 
where 



10 
3 



1 



o(k) 



<52,o(k) 



10 , fe2 

T + IT 



+ 0(r-9/2), (43) 



'52,c(k) = ^ / '^q^i'^(q.k-q) 

x5i,c(q)'5i,c(k- q). 



(44) 



(27r)3 



dq 



3 fc2 



Ukl 



X3i(q)3i(k - q)(5i,c(q)(5i,c(k - q), (45) 



and F^'' is a mathematical function given by eq. [A28| . 
The third-order solution is 



53 (k,' 



7-6 



1 



c(k) 



<53,c(k) 



fc2 

IT 



+ 0(r-"/2) 



(46) 



where 

<53,c(k) : 
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and 



(27r)e 



rfqidq2rfq3(5z)(k - qi - q2 - qs) 



xF3('^(qi,q2,q3)'5l,c(qi)<5i,c(q2)<5i,c(q3), (47) 



dqidq2dq3(5z)(k - qi - q2 - q3) 



(27r)e 



X^i'^(qi,q2,q3)<5l,c(qi)-Jlx(q2)-Jlx(q3), (48) 

(s) (s) 

and and are mathematical functions given by 
eqs. |A30| and [B27| . respectively. One may check that 
these functions arc properly normalized, i.e., gn ^ 1 as 
fc ^ 0, using S'2 c ^ ^i,c and (^3 — ^ (^s.c as /c ^ 0. 
Ignoring the decaying modes, let us rewrite 32 and 33 



as 



52(k) = . 



10 



■52,o(k) 



<52,<:(k) 



-I , 3 fe2 



.93 (k): 



1 



1 fc2 

7 fej 



(49) 



(50) 



These results may be interpreted as, roughly speaking, 
the non-linear filtering functions having smaller effec- 
tive filtering scales (larger filtering wavenumbers): fcj — > 

fcj — \l^kj for the second order, kj ^ kj — \/7kj 



cj = Y 

for the third order, and kj 



k.i - 



|n(n -|- ^)fcj for 

the n-th order perturbations. In other words, the higher- 
order solutions for 5n,b are less suppressed relative to the 
CDM solutions. In the next section we shall quantify 
this effect in more detail. 

3. POWER SPECTRUM 

In this section, we calculate the non-linear matter 
power spectrum using the results obtained in the pre- 
vious section. The total matter fluctuation, 5, is given 
by (5 = fc6c + fb^b, and thus the total matter power spec- 
trum, Ptotik), is given by the sum of three contributions: 

Ptotik^r) ^ f^P,{k,T) + fJMk,T) + f^P,{k,T), (51) 

where Pc{k) and Pb{k) are the power spectra of the CDM 
and another matter component with pressure, respec- 
tively, and Pbc{k) is the cross-correlation power spec- 
trum. Each term is the sum of the linear part, Pii(fc, r), 
and the non-linear parts, P22{k,T) and Pi3(fc,T): 

P^{k,T) = Pll,^{k,T) + P224k,T) +2Pi34k,T), (52) 

where i = (c, 6, be). 

The 3PT power spectrum of CDM has been found in 
the literature (see iBernardeau et ai1l2002l . for a review): 

P22Ak,r)^2 J ^Pii_,(9,r)Fn,c(|k-q|,r) 



Fi^)(q,k-q) 



(53) 



where is a mathematical function given by eq. |A28| , 



dq 



27r 

Pi3Ak,r) = —ePiiAk,r) J ^Pn^c(g,T) 



50 



fc2 

3(g^ 
2 



21 



fc4 

fc2)3(2fc2 





79 



k^q^ 



{2n) 
.k^ 



Ifc- 



(54) 



See Appendix §[X]for the detailed derivations. 

Here, we have implicitly generalized the results from an 
EdS universe to general cosmological models, by writing 



Pii{k,T 



Pll(fc,T, 



DHt) ( 6[%{k,T)/6f:l4k,r) ' 
DHn)[s[%ik,r.)/S^'Uk,r.), 



Pii{k,Ti 



(55) 



where is some arbitrary epoch, r, is the epoch where 
the pressure effect becomes non- negligible (i.e., reion- 
ization epoch for baryons and the relativistic to non- 
relativistic transition epoch for massive neutrinos), and 
D{t) is the linear growth factor appropriate to a given 
cosmological model. This simple generalization has been 
shown to provi de an excellent approx imation to the full 
calculation: see lBernardeau et al.l (1200 21 for models with 
non-zero curvatu re and/or a cosmological constant, and 
iTakahashU ()2008l ) for dynamical dark energy models with 
a constant equation of state of dark energy. 

The linear spectra of the other contributions, Pu^bc 
and Pi 1,6, are given by 



Pll,bc(fc,T) =.gi(fc)Pii,c(fc,T), 
PnAk,r)^gj{k)PnAk,T). 



(56) 
(57) 



The non-linear terms, the main results of this paper, are 
given by 



P22,bc{k,T) 



1 



+ 



10 
3 

14 



[^22,c(fc,T) 



dq 



Piix(g,r)Pii,,(|k-q|,T) 



3 J (2^)3 
xPi^)(q,k-q).Fi^'(q,k-q) 



(58) 



P22,6(fc,T) 



(f +§) 

28 r dci 



i2n) 

l,k 
dq 



2 [^22,c(fc,T) 

3-Pii,c((7,r)Pn,,(|k-q|,r) 



xPi^)(q,k-q).Fi^'(q,k-q) 
98 



27r)3 
(.Fi^)(q,k-q) 



Piix(g,r)Pii,e(|k-q|,r) 
2" 



(59) 
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-Pl3,6c(fc,T)^ 



9i{k) 



A;2 



xPii,c('7,t)] , 



(60) 



.91 (fc) 
7+11 



X / ^^i^^(q,-q,k)Pn,c(g,r) 



(61) 



See Appendix §[C]for the detailed derivations. 

How would Ptot{k) compare with the CDM part, 
Pcik)? 

• In the linear limit, we should recover 
Ptot{k) / Pc{k) -> [fc + fbgi{k)f, which approaches 
unity as fc — )■ 0. 

• In the very small scale limit (fc — > oo), the pres- 
sured component is completely smooth {Sf,{k) 0) 
because gi{k) 0; thus, Ptot{k) / Pc{k) approaches 
a constant value, 

• In the intermediate regime, especially at the tran- 
sition scale between the super- Jeans scale (fc < fcj) 
and the sub-Jeans scale (fc > fcj), the shape of 
Ptot{k) / Pc{k) is significantly distorted away from 
the linear prediction. Non-linear clustering of the 
pressured component adds power at fc '--^ fcj, which 
shifts the effective filtering scale to smaller spatial 
scales as we go to lower redshifts. 

In Fig. [2] we show the ratio, Ptot{k, z)/Pc{k, z) (solid 
lines), for different redshifts (z = 0.1, 1, 3, 5, 10, and 
30), and different fcj (fcj = 1 and 3 h Mpc~^ for the 
left and right panels, respectively). In the linear regime 
(see the bottom lines, z = 30) the ratio agrees with the 
linear prediction shown by the dashed lines. As we go 
to lower redshifts, we find that the filtering wavenum- 
bers continue to shift to larger values, i.e., the filtering 
scales continue to shift to smaller spatial scales as we go 
to lower redshifts. This effect cannot be predicted from 
the linear theory, where all the modes evolve in the same 
way. 

4. COMPARISON WITH APPROXIMATE TREATMENT OF 
SAITO ET AL. (2008) 

The non-linear power spectrum with a significant con- 
tribution from a pressured component has not been 
studie d very muc h in th e literature, with one excep- 
tion. ISaito et all (|2008ll (hereafter, STT) have stud- 
ied effects of massive neutrinos on the no n-linear mat- 
ter power spectrum usi ng 3PT (also see IWond 120081 : 
iLesgourgues et all I2009D . However, their treatment 
is not satisfactory: they have entirely ignored non- 
linearities in neutrinos, but approximated the neutrino 
perturbations as linear perturbations. More precisely, 
they calculated the non-linear matter power spectrum as 

PfJT(fc,z) = f^P,{k,z)+2fJ,PnMk,z)+f.PiiAk,z). 

(62) 




Fig. 2. — Ratio of the total matter power spectrum, Ptot{k,z), 
to the CDM part, Pc(fc,z), at z = 0.1 (top), 1, 3, 5, 10, and 30 
(bottom). (Left) The input Jeans wavenumber of fcj = 1 h Mpc^^. 
(Right) kj = 3 h Mpc~^. The dashed Unes show the ratios calcu- 
lated from the linear theory, whereas the dot-dashed lines show the 
linear calculations with fcj = 2 and 6 h Mpc~^ for the left and right 
panels respectively, to show that the actual filtering wavenumbers, 
predicted by the 3PT calculations, can be ~40% as large as the 
linear filtering wavenumber at low redshift. 

In our language this leads to 



jSTT 



{k,z)^f,P4k) + [2fJ,gi{k)~ 



f,9l{k)]PnAk,z). 

(63) 

■'b" with 'V" to 



Here, we have replaced the subscripts 
avoid confusion in notation. 

How accurate is the STT approximation? To study 
this, we compare eq. |63| to the full calculation given in 
the previous section. Figure [3] shows the fractional dif- 
ference between our full calculation and STT's approxi- 
mation, [Ptot{k) - PgY (k)]/ Ptot{k), for njn^ = 1/10, 
1/20, and 1/100, which correspond to the sum of neu- 
trino masses of J2i — 1-3, 0.64, and 0.13 eV, respec- 
tively, where i = {e,iJ,,T). We find that STT's approx- 
imation clearly underestimates the power at fc ~ Afs, 
where kps is the neutrino free-streaming scale, or it is 
the Jeans wavenumber computed with the velocity dis- 
persion of the neutrinos. More precisely. 



kps.ii'T) 



V6 



(64) 



in an EdS universe, where cr^j(T) is the velocity dis- 
persion of neutrino species i (see, e.g.. Appendix A. 3 of 
iTakada et aHl2006D . 

One may argue that STT's approximation should be 
better for a smaller neutrino mass: the errors in the to- 
tal matter power spectrum are 3.5%, 0.6%, and 0.003% 
X^i "^i^.i ~ l-^j 0.64, and 0.13 eV, respectively, at z = 0.1; 
however, our results indicate that their approximation is 
conceptually not correct: neutrinos should not be treated 
as linear perturbations, as the neutrino velocity disper- 
sion has no effect in suppressing the neutrino perturba- 
tions at and above the free-streaming scale. In other 
words, the errors may happen to be small in the total 
matter power spectrum for small neutrino masses be- 
cause neutrinos contribute only a tiny fraction of the to- 
tal matter density anyway, but the errors in the neutrino 
power spectrum are large. Figure |4] shows the fractional 
difference between the non-linear neutrino power spec- 
trum, P^(fc), and the linear power spectrum, P^*"(fc), 
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Fig. 3. — Fractional diffe rence between our full calculation and 
the approximation used by ISaito et"all ll2008i) (STT), [Ptot(k) - 
PfJt'^(k)]/Ptot{k), for C^/n^ = 1/100 (top), 1/20 (middle), and 
l/lO (bottom), which corresponds to '^rrii, ~ 0.13, 0.64, and 
1.3 eV, respectively. 

i.e., AP/P = [Pu{k) - Pl''^{k)]/P^{k). It is clear that 
neutrinos are significantly non-linear, even well below 
the free-streaming scale, k ^ kps- Nevertheless, the 
STT approximation may still provide a convenient phe- 
nomenological tool for calculating the non-linear total 
matter power spectrum in the presence of massive neu- 
trinos. 

5. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have obtained the second- and third- 
order solutions for the density perturbations in a system 
consisting of two matter components with and without 
the pressure gradient. This is the first self-consistent an- 
alytical calculation, with non-linearities in the pressured 
component fully retained up to the 3rd order in pertur- 
bations. 

As our study is focused on understanding the physics of 
the non-linear pressure effect on the matter power spec- 
trum, we have studied a toy model in which the Jeans 
wavenumber, fcj, is independent of time. This is equiva- 
lent to the temperature of the pressured component fol- 
lowing that of radiation, i.e., T oc 1/a. 

Nevertheless, we have found several results that have 
qualitative implications for the practical applications. 
We have found that non-linearities in the pressured com- 
ponent shift the f iltering scale from th e well-known lin- 
ear filtering scale (|Gnedin fc Hui|[l998[ ) to a smaller spa- 
tial scale (larger wavenumber) by a factor depending on 
the redshift and the Jeans scale. In other words, the 
actual filtering scale for a given sound speed (or tem- 
perature) is smaller than the linear scale. Therefore, if 



Fig. 4. — Fractional difference between the non-linear neutrino 
power spectrum, P^(fc), and the linear power spectrum, Pjf"{k), 
[Pu(k) - Pl''\k)]/ pS{k), for n^/n™ = 1/100 (top), 1/20 (middle), 
and 1/10 (bottom), which corresponds to '^rriu ~ 0.13, 0.64, and 
1.3 eV, respectively. 

one used the linear filtering scale to interpret the fall-off 
of, e.g., the flux power s pectrum of the Lyman-a forests 
(jZaldarriaga et aLl l2001h. one would underestimate the 
temperature of the pressured component. 

How important is this effect? For example, when the 
Jeans wavenumber is kj — IQ h Mpc~^, our calcula- 
tion predicts that the effective filtering wavenumber is 
~ 10, 12, 13, 13, and 14 h Mpc'^ at z = 30, 10, 5, 
3, and 1, respectively. While we do not expect 3PT to 
be valid at such high wavenumbers, our results clearly 
indicate that the expected changes in the filtering scale 
cannot be ignored. Table [T] summarizes the ratios of the 
effective (actual) and the linear filtering wavenumbers. 
Note that the linear filtering wavenumber is the same 
as the Jeans wavenumber in our model; thus, we show 
kp^eff/kj in Table[TJ We extracted the effective filtering 
wavenumber, kF,eff, by fitting [fc + fb/{l + k"^ / kjr^ef f)? 
to Ptot(k, z)/Pc{k, z). We find that a factor of 1.4 change 
in the filtering scale is quite common over a wide range 
of redshifts and kj. 

A factor of 1.4 change in the filtering scale changes 
the inferred temperature by a factor of two; thus, one 
implication of our result is that the temperature of the 
Inter-galactic Mediu m (IGM) obtained from the Lyman- 
a forests at z = 3 bv lZaldarriaga et all ()200lD might have 
been underestimated by a factor of two. 

A factor of 1.4 change in the filtering scale gives a fac- 
tor of ^ 3 change in the filtering mass. Our calculation 
shows that the actual filtering mass is similar to the lin- 
ear one only in high redshifts, while the former is signifi- 
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TABLE 1 

Ratio of the effective and the linear filtering 

SCALES, kp^eff/kj 





z=0.1 


1.0 


3.0 


5.0 


10 


30 


[h Mpc'i) 














0.1 


1.08 


1.04 


1.01 


1.00 


1.00 


1.00 


0.5 


1.37 


1.21 


1.07 


1.03 


1.01 


1.00 


1.0 


1.43 


1.32 


1.14 


1.08 


1.03 


1.00 


3.0 


1.41 


1.38 


1.28 


1.20 


1.08 


1.01 


5.0 


1.40 


1.39 


1.32 


1.24 


1.12 


1.02 


10 


1.41 


1.40 


1.35 


1.29 


1.16 


1.03 



Note. — This table shows the ratios of the effective 

//) 
lifte 



(fc^ gjf) and the linear (fcj) filtering scales for different 
redshifts and kj. The ratios are closer to unity at higher 
redshifts because non-linearities are weaker. 

cantly smaller than the latter in low reds hift. This result 
is qua litati vely similar to thos e found in lOkamoto et all 
(|2008D and iHoeft et all (|2006[) : however, a quantitative 



comparison is not possible, as our results apply only to 
the system with a constant Jeans wavenumber. 

What is next? As for baryons, we need to extend our 
formalism for incorporating a realistic thermal history of 
the universe with a proper time dependence of kj. As 
for neutrinos, we need to incorporate not only the pres- 
sure gradient but also the anisotropic stress in the Euler 
equation. To do this we need to solve the Boltzmann 
equation. Nevertheless, our results presented in this pa- 
per show that neutrinos are significantly non-linear, even 
well below the free-streaming scale. 

This material is based in part upon work supported by 
the Texas Advanced Research Program under Grant No. 
003658-0005-2006, by NASA grants NNX08AM29G and 
NNX08AL43G, and by NSF grant AST-0807649. E. K. 
acknowledges support from an Alfred P. Sloan Research 
Fellowship. 



APPENDIX 
3PT FOR CDM 

The continuity, Euler, and Poisson equations of CDM are given by 
• Continuity equation: 



• Euler equations: 



<5(x,T) + V-[(l + <5(x,T))v(x,r)] -0, 



v(x, t) + [v(x, t) • V]v(x, r) = ^ '^<^(^' 



(Al) 
(A2) 

(A3) 

First, we take the divergence of Eq. |A2| and substitute Eq. |A3j . Moving all the non-linear terms to the RHS of the 
equations, we find 



Poisson equation (for an EdS universe): 



V2</,(x,r) = -<5(x,r). 



^(x,r) + V • v(x,t) ==-V- [(5(x,T)v(x,r)], 

Qi 6 

— [V • v(x, r)] + - [V • v(x, r)] + — (5(x, r) = -V • {[v(x, r) • V]v(x, r)}. 
OT a 



Let us take the Fourier transform of Eqs. |A4| and |A5I 

1 



J(k, r) -f 0(k, r): 



(27r)3 



dqidq2(5D(qi + q2 - k) — —e{qx,T)5{ci2,T), 

9i 



(A4) 
(A5) 

(A6) 



^^(k,T) 



1 



dqidq2(5D(qi + q2 - K 



fc^(qi ■ q2)2 



a - (27r)3 

where we have defined ^ = V • v, and its Fourier transform is given by 

k ~, 



^^(ql,r)0(q2,r), 



v(k,r) = -*^^>,T). 



(A7) 



(A8) 



One can decompose the solutions of the non-linear continuity and Euler equations, S and 9, into the sum of infinite 
series of n-th order perturbations of density and velocity divergence fields: 



^(k,r) = ^a"(r)<5„(k), 

oo 

0(k,r):=^<i(T)a"-i(T)0„(k), 



(A9) 



(AlO) 



respectively. Note that, strictly speaking, this particular decomposition, a decomposition into a series with powers of 
a(T), is valid only for an EdS universe. However, generalization to arbitrary cosmological models can be done in the 
end by replacing air) with the appropriate growth factor, D{t) (|Bernardeau et aLll2002l : lTakahashill2008f ). 
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by 



Now, let us solve Eqs. |A6| and |A7| at each order of perturbations. The n-th (n > 1) term of the Eq. |A6| is given 



d(T)a"-i(T)H„(k)+0„(k)] =- 



(27r)3 



dqidq2(5n(qi + q2 -k) ^ ^ d(r)a" ^(T)6'„i(qi)(5„_m(q2)(All) 

9i ^=1 



n-l 



Dividing both sides by d(r)a" ^(t), one obtains 



where 



(27r) 



Similarly, from the Euler equation, Eq. |A7j . one obtains 

35„(k) + (l + 2n)f?„(k) -B„(k), 

where 



Bn(k) 



1 



(27r) 



n(5„(k)+0„(k) = ^„(k). 



dqidq2(5z)(qi + q2 - k) — j-!- ^ 6',„(qi)(5 

11— rn (q2). 

^1 m=l 



ri-1 



dqidq2dD qi + q2 - k) — 

9r'?2 



(q2). 



(A12) 
(A13) 
(A14) 
(A15) 



m— 1 



The forms of Eqs. [A12| and [A14| indicate that the n-th order solutions are written in terms of the sum of 1-st to 



(n — l)-th order solutions, with Si(k) 



^^i(k). By solving Eqs. |A12j and |A14j for (5„ and 9n, one obtains 
(l + 2n)A„(k)-B„(k) 



<5n(k) = 
0„(k) 



(2n + 3)(n- 1) 
-3A„(k) +ng„(k) 
(2n-|-3)(n- 1) ' 



which can be rewritten as 
<5n(k) = 

e„(k) = - 



1 



(27r)3"-3 
1 



dqi...dq„^D(qi + ••• + q« - k)F„(qi, ...,q„)(5i( qi j.. .01 (qn j , 
dqi...dq„5D(qi + ••• + q« - k)G„(qi, q„)(5i(qi)...(5i(q„). 



(27r)3"-3 

Here, the newly defined kernels, F„ and G„, can be found from the following recursion relations: 



n-l 



p / \ _ G'm(qi, qm) 

-f'rUqij qnj — / , /D I o^/' l\ 

(2n -f 3)(n - 1) 



n , o ^1 f N^^Cqi • q2)^ 

U + in) 5 — l'n-m\f\m+l, q^j 2~2 



(,qm+i 



1 • ■ ■ I j 



and 



n-l 



G„(qi, ...,q„) = ^ 



Gm(qi, qm) 
(2n-h3)(n- 1) 



,k - qi 
' T~ 



-^n — m (Qm+1 ; ■■■jQn ) + n 



fc^(qi ■q2) 

? 2 
9192 



Gn — m(qrfj+l, ...,q7i) 



with the boundary conditions of -Fi = 1 Gi. The 2nd-order solutions are 



^2(qi,q2)==3 



5k-qi k'^(s\x ■ ^2) 



7 9f 



+ 



^ , , 3 k • qi 2fc^(qi -qa) 
G2(qi,q2)==7; — — 



7 



7g?gi 



where k = qi + q2 . The 3rd-order solutions are 



■F3(qi,q2,q3) 



1 

18 



-^^^^^2(q2, qs) + ^-^7572^^G2(q2, qs) 



11 



2 2 
9l923 



G'2(qi,q2) 
18 



7k • qi2 fc^(qi2 • qa) 



9?2 



2 2 

91293 



where q^ j = q^ -t- q^ and k = X] qi ■ 
It is often convenient to have the symmetrized forms of the above kernels. They are 

-F'i'''(qi>q2) = ^ [^2(qi,q2) + ^2(q2,qi)] , 
^^'■'(qi, q2) = ^ [G2(qi, q2) + G2(q2, qi)] , 

(qi , q2 , qs) = ^ [^3 (qi , q2 , qa) + ^3 (qi , q3 , q2) + (pa , qi , q3) 
+ -F3(q2,q3,qi) + -F3(q3,qi,q2) + ^3(q3,q2,qi)] 



(AI6) 

(A17) 

(A18) 
(A19) 

, (A20) 

, (A21) 

(A22) 
(A23) 

, (A24) 

(A25) 
(A26) 

(A27) 
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The explicit forms are 



G^''(qi,q2)^ 



i^i'^(qi,q2,q3) = -k- 



5 2 (qi.q2)" l (qi-q2)(gHgi) 

7 7 +2 

3 4 (qi ■ q^)^ , 1 (qi • <12M + ql) 

7^7 qlql +2 
7 

54 



9 ? 



(A28) 
(A29) 



f!,'^ (q2 , qs) ^ + Fi^) (qi , q3) ^ + F^^^ (qi , q^ ) 3| 
9i I2 % . 



r^[s) , Nqi-q23 , r^{s), Nq2-qi3 , ^(s), n 
^2 (q2,q3) ^2„2 +^^2 (qi,q3) , , +G'^ ^(qi,q2)— ^-2- 



, q3 • qi2 



91923 



9593-3 



^(s), ^q23 , ^,{s), sqi3 , ^(s)/ \qi2 
G2 (q2,q3)— + (qi,q3)— + (qi,q2)— 

923 9l3 9l2 , 



(A30) 



In order to calculate the next-to-linear-order density power spectrum, one needs to use the solutions of the density 
fluctuations up to the 3rd order: 

(2^)3p(fc, t)6d{\^ + k') = (5(k, T)^(k', r)) 



kI=1 



which yields 



.a2(r)(5i(k)<5i(k')) +a^(r)(<5i(k)53(k') +'52(k)<52(k') +<53(k)<5i(k')), 



P{k, t) ^ a\T)Pn{k) + a^T) [P22(fc) + 2Pi3(A:)] + OiS"^). 



Here, we have defined the quantity, Pij{k), given by 



The non- linear corrections, P22{k) and Pi^{k), are 

dq 



P22(fc)-2 



Pn(9)Ai(|k-q|) Pi^)(q,k-q) 



(27r)= 



where 



and /i = k • q, and 



^^i'^(q,k-q) 



7 14 



q'^{k'^ 



2kqii) 



Pi3(fc) = 3Pn(fc) 



^\pi^)(q,-q,k)Pn(9). 



(27r) 



Using 



^d/xPi^)(q,-q,k) = — 



5 - 2 1 - - 79 S + 6 ^ + ? i^l^^Sl±M 1, , ^ + 



9 



k^q^ 



\k-q\ 



one obtains (jMakino et al.|[l993 ) 



soil - 21I; - 79 + 6^, + ^ in ^ + ' 



fc2 



|fc-9l 
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(A31) 
(A32) 

(A33) 

(A34) 

(A35) 
(A36) 



(A37) 



In this Appendix we shall derive the higher-order filtering functions. We shall solve Eqs. |14]-[17| perturbatively, up 
to the 3rd-order in perturbations. The density contrasts and velocity divergence fields of CDM and the matter with 
pressure are all expanded into the infinite sum of n-th order perturbations as 



5e(k,r) = ^a"(r)5„,,(k), 

00 

^e(k,T) =^d(T)a"-l(T)0„,,(k), 



(Bl) 
(B2) 
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5b(k,r) = ^a"(r)5„,,(k)5„(k,r), 

oo 



(B3) 
(B4) 



where g„(k, r) and /i„(k, t) are the filtering functions for the density and velocity divergence fields, respectively, at 
the n-th order. 
With the above series expansion, Eqs. [15] and [17] yield 



[(na(r)a"-i(r)g„(k,r) +a"(T)5„(k,T)) <5„,,(k) + a(T)a"-i (r)/i„(k, r)0,,,(k)] 



(27r)3 



dqidq2(5D (qi + q2 - k) — ^ ^ ^ i 

'^1 m,=l ;=i 



^^m(qi,T)5;(q2,T)6'mx(qi)<5/x(q2), 



(B5) 



oo 

^ [(a(r)a"-i(T) + a2(T)a"-2(r)(n - 1)) /i„(k, r)0„,,(k) + a(r)a"-i(r)/i„(k, r)0„,,(k) 

ri=l 

+ -a(T)a"-i(T)/^„(k)(k, T)0„,e(k) + (/c + A5n(k, r)) ,5„,,(k) 

A n p 1,2 f xoooo 

= / / rfqirfq2'Jz.(qi +q2 -k)^^^^ E E'^'W«'"'''"'M'^™(qi'^)'^'(q2,r)0„,,(qi)0i,,(q2) 

OC 

+fc2c2(r)^a"(T)g„(k,T)^„,,(k) 

x3 fe^cg(T) / / dqidq2(5i5(qi +q2 -k) ^ ^a'"+'(T)5„,(qi,r).g,(q2,T)(5„,c(qi)(5/,c(q2) 

^ ^ m—1 l—l 



dqi(iq2dq3(5n(qi + q2 + qs - k) 



OO OO Ctt 



^ 5Z XI XI a'"^'^^W5m(qi, r).g;(q2, T).gp(q3, r)(5™,c(qi)5ix(q2)5p,c(q3)- 

m—1 l — l p—1 



(B6) 



From now on, we shall write the sound speed, Cs, in terms of the usual Jeans wavenumber, kj, as Cs = VQ/ {kjr). We 
shall ignore the inhomogeneity in Cs_ (i.e., spatial dependence of Cg) throughout this paper. For the linear analysis for 
Vcs 7^ 0, see lNaoz fc Barkan 



Second Order Solutions 

We have derived the linear filtering function, gi{k), in Eq. [30j . For n = 2, the continuity and Euler equations are 
given by 

<52.c(k).92(k,r) + -52,c(k)52(k, r) + -02,,{k)h2{k,T) 

T T ' 



2 1 



T (27r)3 
-^2(k) 



dqidq25Li(qi + q2 - k) — 2~'^i.c(qi)'^i.c(q2)5i(qi)ffi(q2) 

9i 



(B7) 



10 2 6 6 fc^ 

— 02,c(k)/l2(k,T) + -02,c(k)/l2(k,T) + — <52,c(k) - ^^<52,c(k)g2(k,T) 



r2 ^2 



4 1 
72 (27r)3 

^i?2(k) 



dqidq2(5D(qi + q2 - k) 



3fc2 fc2(qi.q2) 



4 fc^ 2q2g2 



'5i,c(qi)^ix(q2)5i(qi).9i(q2) 



(B8) 



Here, 6*1, c(k) = — (5i,c(k). Combining Eqs. [B7| and |B8j . we get the second order inhomogeneous partial differential 
equation: 
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•• /, ^ 10 . , 1 

92 (k, t) + —g2 (k, t) + — 



20 + 6 



ff2(k, r) 



10^2 (k) , 4^2 (k) 



<52,c(k) <52,c(k) 



(B9) 



where (52,c(k) is given by 



^2,c(k) = ^ / ^q^i'^(q>k-q)(5i,,(q)Ji,,(k-q). 



(BIO) 



Solving the above differential equation, we have: 



10A2(k) 4B2(k) 



20 + 6 



(Bll) 



where the oscillation component, 



1±. 1 



(B12) 



decays for any choice of < k/kj. The 2nd-order filtering function for the velocity divergence field, /i2(h, r), is given 

by 



h2{k) = ^-Ij^ [^2(k) - 2(52,c(k)52(k)] , 



(B13) 



where we have ignored the decaying term. 
Using the explicit forms of A2(k) and -B2(k) given by Eqs. |B7j and |B8j . respectively, we obtain 



52 (k, r) 



10 _ 7 
3 3 



■52.c(k) 



10 I fc2 

T + IT 



+ 0(r-9/2), 



(B14) 



where (54 „ is 



'^2,c(k) = ^ / ^q-^i'^(q'k-qMi,c(q)<5i,c(k-q), 



(B15) 



where 



-^2 (qi.q2) = 



^i^^(qi,q2) . ^4^2 



3 fc2 



5i(qi)ffi(q2)- 



(B16) 



In the limit where kj oo, T2'^\qi, q2) = F2''''(qi, q2), and thus 172 ^ 1- For the velocity divergence filtering function, 
we find 



/i2(k)= 



1 



5^2,c(k) 
1 



72, c 



(k) 



~^JJ dq_idq2SD{q.i + q2 - k)(5i,c(qi)5i,c(q2)(2F2'''-'(qi,q2)- G'^'''(qi,q2)).gi(qi).gi(q2)-252x(k)52(k) 

'52,c(k) 



(27r) 
1 



i2ny 



dqidq2(5i5(qi + q2 - k)(5i,c(qi)(5i,c(q2) ( H — f^-^ — )5i(qi)ffi(q2) 



-2 



^2,c(k) 



52(1^1817) 



where we have used 2F2(qi, q2) — G2(qi, q2) = This expression also converges to /i2 = 1 as we take the limit of 

kj — > 00. 
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Third Order Solutions 
For n = 3, the continuity and Euler equations are given by 

3d(T)a2(T)53(k, r),53^c(k) + a\T)g^{]^, r)<53.c(k) + h{T)a\T)h^{\,, T)03,c(k) 



1 



(27r)e 



dqidq2dq3(5i3(qi + q2 + qs - k)(5i,c(qi)(5i,c(q2)^i,c(q3) 



^ ^^5i(qi)ff2(q23)-Fi''''(q2,q3) + ^^^^^2(qi2)5i(q3)G^'*Hqi>q2) 



ll2 



a(T)a2(r)A3(k), 



[S(r)a2(r) +2a2(r)a(r)] /i3(k, r)03,c(k) + a(T)a2(r)A3(k, T)03,c(k) + -a(T)a2(r)/i3(k, r)03,c(k) 
+ 4«'W'53,c(k) - 4Sa'(r)<53,c(k) 



r2 



(27r) 



' ^-T^TT^ff 1 (qi ) ^2 (q23 ) G^"*^ (q2 , q3 ) 

^91 923 



dqidq2dq3<^n(qi + q2 + q3 - k)(5i,c(qi)(5i,c(q2)^i,c(q3) 

fc^(qi2 ■ q3) 



29?29i 



-/i2(qi2)gi(q3)G^'*^(qi,q2) 



3 fc^ f -1 3 A:^ f ^ 1 A:^ 

4p-5i(qi).92(q23)-F2 (q2,q3) - |p-.g2(qi2)5i(q3)^2 (qi.q2) + 2p"5i(qi)5'i(q2)5i(q3) 



_ -2 



(r)a(T)i?3(k). 



In an EdS universe, a^r) = we have 



^3.c(k)g3(k,T) + -<53.c(k)g3(k,T) + -03.c(k)/i3(k, r) - -A3(k), 

T T T 

^/i3(k,T)03,c(k) + -A3(k,T)03,c(k) + 4'53,c(k) - -|S-'53,c(k)g3(k,T) = 4s3(k). 

7" 7" T T rv J T 



Combining Eqs. |B20j and |B21| . we have the second-order differential equation: 



14 1 
53(k,T) + — 53(k,r) + — 



42 + 6-^ 53(k,r) 



14^3 (k) 4i?3(k) 

<53^c(k) ^ (53^c(k) 



= 0. 



Solving this, we obtain 



^ 7A3(k) _ 2B3(k) 
g3(k, r) = 3.3.(k) 3.3.(k) ^ 0(^-13/2)^ 

7+7:5- 



(I 



(I 

(I 
(I 

(I 

(I 



-13/2(l±./l-^(7-Hf|-) 



where the oscillation component, 

0(r-i3/2) 0, 

decays for any < pj- . The velocity divergence filtering function at the 3rd-order is 

^3(k) = tAt [^3(k) - 3<53,c(k)53(k)] , 



^3,c(k) 



where we have ignored the decaying term. 
Let us rewrite 7^3 (k) - 2B3(k) in Eq. [B23j as 



7A3(k)-2B3(k): 



1 



(27r)6 
"7k qi 



9i 



dqidq2dq3(5£,(qi + q2 + q3 - k)5i,c(qi)(5i,c(q2)'5i,c(q3) 



.9i(qi)ff2(q23)-Fi'-'(q2,q3) + -^^^-^ft-2(qi2)gi(q3)G^'*^(qi,q2) 



9?2 



+ 



3i(qi)/i2(q23)G^'*^(q2,q3) 



fc^(qi2 • q3) 



? 9 

91293 



/i2(qi2).gi(q3)G^'*^(qi,q2) 



fc^(qi ■ q23) 

2 ? 

91923 



(I 
(I 



+ o 712-51 (qi)32(q23)^2 (q2,q3) + 2p-£'2(qi2)5i(q3)i^2 (qi.q2) - i:2 5i(qi)ffi(q2)5i(q3) 



2A:2 
18 

18^3,c(k) 



dqidq2C?q3(5n(qi + q2 + q3 - k)J^3(qi, q2, q3)(5i,c(qi)<5i,c(q2)(^ix(q3) 



(I 
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The new kernel, .F3(qi, q2, qs), can be symmetrized as 

-^3'*^(qi,q2,q3) 

= ^ [^3(qi,q2,q3) +^3(qi,q3,q2) +^3(q2,qi,q3) +^3(q2,q3,qi) +^3(q3,qi,q2) +^3(q3,q2,qi) 



■f^i'''(q2,q3)%gi(qi)52(q23) + i^i"-'(qi,q3)-|5i(q2)32(qi3) + i^i'-'(qi,q2)%ffi(q3)52(^^ 
ii I2 I3 

G'2''^(q2,q3) % 2^^ ffl(qi)fe2(q23) + G'2''^(qi,q3) ^'^2 2" .9i(q2)^2(qi3) + <?2^^(qi,q2) ^\ ^^^ gi(q3)/t2(qi2) 

9l923 92 913 93 912 



7 



k- 



+ 



54 
18 kj 



G^'^(q2,q3)^ffi(qi)/i2(q23) + G^''^(qi,q3)^5i(q2)/i2(qi3) + G^'^(qi,q2)-'^5i(q3)/i2(qi2) 

923 9l3 9l2 



5i(qi)52(q23)i^2^\q2,q3) +5i(q2)52(qi3)i^2''^(qi,q3) + 5i(q3)52(qi2)i^2^^(qi,q2) -5i(qi)5^ 



(B27) 



In the limit of kj ^ 00, ^ F3, and gsik) = 1. Using ^(k) introduced above, we write (?3 as 



53(k) = 



7-6 



1 



7 + 



fe2 



(B28) 



3PT TOTAL POWER SPECTRUM 

We calculate the power spectrum of the total matter fluctuations, 5 = fcSc + fbh = fcSc + (1 — fc)Sb, which is given, 
up to the third-order in perturbations, by 

6(k,T)=fMKT) + fbSbiKr) 

= [(5i,e(k, r) + 52,c(k, r) + (53,c(k, r)] + (1 - /,) [<5i,6(k, r) + <52,b(k, r) + 53,b(k, r)] 

= /e [5i,c(k, r) + ,52,c(k, r) + Ss,c{^, r)] + (1 - A) [ffi(A;)5i,e(k, r) + ff2(k)<52,c(k, r) + 93(^)63 r)] . (CI) 
The power spectrum is 

{2TrfPtot{k)dD{k + k') = ((5(k)(5(k')) 



= ({/A(k) + (1 - /c)56(k)}{/e^c(k') + (1 - /c)^6(k')}) 

= /2(<5,(k)<5,(k')) + 2/e(l - /e)(56(k)5e(k')) + (1 - fcf {Sb{m{^')) 

= (27r)3 [/2p,(fc) + 2/e(l - fc)PbAk) + (1 - fcfPbik)] 5zj(k + k'), 



(C2) 



where Pc, Pb,c and Pb are 
{27rfPc{k)SD{k + k') = (<5e(k)5c(k')> 



: ({<5i,e(k) + d2,c(k) + 53,c(k)}{5i,e(k') + (52,c(k') + dsA^')}) 
:(5l,e(k)5i,e(k')) +2(5i,e(k)53,c(k')) + (52,c(k)52,c(k')) 

: (27r)3 [Pii,c(fc) + 2Pi3,c(fc) + P22,c(fc)] SdO^ + k'), 



(C3) 



{2TrfPbc{k)5D0^ + k') = {5b(k)S,(k')) 



-- ({5i(fc)(5i,e(k) + ,g2(k)<52x(k) + g3(k)<53x(k)}{5i,,(k') + (52,c(k') + S^A^')}) 

= 5i(fc)(<5i,c(k)(5i,e(k')> + ,gi(fc)(<5i,c(k)(53,e(k')> + (g3(k)(53,c(k)5i,c(k')> + (<?2(k)(52,c(k)52,c(k')) 

I {2TTf [Pn,bc{k) + 2Pi3,bc(fc) + P22,fac(fc)] Soik + k'), (C4) 



{27rfPbik)6D{k + k') = {Sb{k)Sb{k')) 



= ({ffl(fc)5i,e(k) +ff2(k)52,c(k) +53(k)53,c(k)} X {5i(fc')5l,c(k') +52(k')^2,c(k') +53(k')^3,c(k')}) 

■.gl{k){S,Ak)5iAk')) + 2.9i(fc)(5i,,(k)g3(k')<53,c(k')) + (52(k)52,c(k)52(k')52,c(k')) 

; (27r)3 [Pii,6(A:) + 2Pi3,fe(A:) + P22,6(A:)] fe(k + k'), (C5) 



respectively. 
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Wovenumber [h Mpc"'] 
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Fig. 5. — The dimensionless power spectra, A^(fc) = k^P{k)/{2n^), for a matter component with pressure (i.e., baryon, neutrino, etc) 
are shown for several redshifts {z = 0.1, 1.0, 3.0, 5.0, 10 and 30). We show the non-Unear calculations with 3PT in the solid and dotted 
lines for kj = 1.0 and 3.0 h Mpc~^, respectively. We also show the linear calculations in the dashed and dot-dashed lines for kj = 1.0 and 
3.0 h Mpc~^, respectively. 
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SHOJI & KOMATSU 



Now, Pii^c(fc), Pi3,c{k) and P22.c{k) can be numerically calculated with the corresponding kernels, i^2'''' and Fg*-*; 



{27rfPiiMk)6Dik + k') = (5i,b(k)Ji,,(k')> 

= 5i(fc)(,5i,,(k)Ji,,(k')), 

(27r)3pi3,,,(A)(5^(k + k') = i [((5i,b(k)(53,c(k')> + (5i,c(k)53,b(k')>] 

= ^ [5i(fc)(^i,c(k)(53,.(k')) + (<5i,,(k)53(k')^3,c(k'))] 



(C6) 



fe2 



(27r)3p22,6c(fc)^D(k + k') = (<52,b(k)J2,c(k')) 

= (52(k)<52,,(k)J2,c(k')> 

1 



10 I 

T + feT 



(5l,e(k)<53.c(k')) + — ^(5l.c(k)(5' (k')) 



(C7) 



(<52.c(k)<52,c(k')) + -{SUk)S2A^')) 



i2TTfPii,t{k)6D{k + k') = (5i,b(k)^i,f,(k')> 

= .9?(fc)(5i,,(k)<5i,,(k')), 

(2^)3Pi3_,(fc),5c(k + k') = (<5i,b(k)53,6(k')) 

= 5i(fc)('5i,c(k)g3(k')53,c(k')) 
,91 (fc) 



fe2 



[(<5i,,(k)<53,c(k'))+6(^l,e(k)5^,,(k'))] 



(C8) 
(C9) 

(CIO) 



(2^)3p22,6(fc)'5l5(k + k') = (<52,b(k)52,6(k')) 



:(.g2(k)(52,c(k).g2(k')^2,c(k')) 

1 



10 I 



(^2,c(k)<52,c(k')> + y(<52,c(k)5L(k')> + f (<5L (k)^L (k')> 



(Cll) 



The ensemble averages of the products involving ^(k) are given by 

(5i,e(k)(5^,,(k'))=3^z5(k + k')Pii.c(fc) ydq^i^Hq,-q,k)Pii.c(g) 

= 67rSD{k + k')PnAk) dq q^PnAl) J ^ ^M-^i'^ (q, -q, k), 

{S^A^Kci^')) = 2SD{k + k') J dqPn,c(g)Pllx(|k - q|)P<^^ (q, k - q)^^'^^ (q, k - q) 



{d',A^)S',A^'))=2dD{k + k') J dqPn.c(g)Pii,c(|k-q|) [^i^)(q,k-q) 
Here, the term, / ^^J^3(q, -q, k)Pii_c(g), in Eq. |C12| is given by 
/ (^•^i'^(q'-q'k)Pii,c(g) 



(C12) 
(C13) 
(C14) 



(2^)3 
dq [7^, 



(27r)3 [54 



P2^^^ (-q, k) ^51 (q)52 (k - q) - ^^2^^^ (q, k) ^51 (q)ff2 (k + q) 



27 

14, 
H k- 

54 



^^^^ ?(k^-^^ (^^■'^^ " " ^^^'^ ?(kT#^^ ^^^^^ + 

F(^) (-q, k) -ii^.gi (q).g2 (k - q) + P*^^ (q, k) ^^^91 (q)52 (k + q) 



, q ■ (k + q) 



-fc2 

7 

-k- 



1 

'27 

7 

"54 



'(k-q)2 
(-q.k)(g2 + fc2)^ q ■ (k - q) 
g2(k-q)2 



2g2fc2 

-,2 I I,2> 



.gf(q)5i(k)- 1 + 



(k + q)2- 

(q.k)(g2 + P)\ q.(k + q) 



2g2fc2 



(k - q)2 



2g2A:2 
2 I ^,2^ 



g2(k + q)2 



5?(q)5i(k) 



^ , (-q-k)(g2 + fc2) ^ k-q 2, ^ ^ /^i ^ ISJlK?_t^^ Ji±3_ 2, , 
1 + I 71 ^5i(q)5i(k) + I 1 + I — — ^5i(q)5i(k) 



2g2/c2 



(k + q)2 



gi(q).g2(k - q)Fi''(-q, k) + gi(q).g2(k + q)F2^''^ (q, k) - .g^(q)5l(k) } Pii,c(g), 



is). 



(C15) 
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where we have used Eq. [B17| and F2*''(q, — q) = 6*2" ''(q, — q) = 0. We then calculate the angular average of JF, 
/ dfj,J-^^\ for the linear filtering function of gi{k) — 1/(1 + k^/kj): 
»i 



3 ' 



pOr^s-^ [-140005" + 810r^"(l + s 



„i 612360r8s(l + r2)(r2 +s , 

60r''s2(i05 - 125s2 + 78s^) + 9r*(321 - 2485^ + ISOs"*) + 27r"(126 - 87s^ + 70s 

1 + s 



900r2s4(-7 + 5s2) 
9,s6)] 



-243r«(-7 + 5s^ + 2s^)[5(r* + s^)(-l + s'Y + r^(5 - 5s^ - 19s^ + 5s'')] In ■ 



1 



2\2 



+ [10s2 + 3r2(l + s2)][-35s^ + ?,r\-7 + s2)][-2000s'' + 135r**(-l + s^) 
+240r^s2(3 - 4s2 + 3s^) + 300r2(s* + s") + 27r'^(5 + bs^ - 9s* + 55*^)]^ In 



10s2 + 3r2(l + s)^ 
10s2 + 3r2(l - s)2 



(C16) 
0). The 



where r = k/kj and s = fc/^/. We find that the calculation of is numerically unstable as k/kj — > (r 
exact limit of is limj,/j,^_^o Tj, F3, and thus one may replace J^^ with F3 for a sufficiently small value of k/kj. 
Finally, we generalize the above results from an EdS universe to general cosmological models, by writing 



Pll(fc,T,)^Fll(fc,T) = 



D'{r) ( 5[%{k,r)/5^^l{k,T) 



A%ik,r^)/S^^}4k,r.)^ 



Pii{k,Ti), 



(C17) 



where is some arbitrary epoch, is the epoch where the pressure effect becomes non- negligible (i.e., reionization 
epoch for baryons and non-relativistic transition for massive neutrinos), and D{t ) is the linear growth factor appropriate 
to a given cosmological model. We obtain Eq. [55] from combining Eqs. |C8| . |C13j . and P22.C given by Eq. A34 . 
Similarly, we obtain E q. [59] from combining Eqs. [Cllj . [C14| . and P 22.c, E q. [50] from combining Eqs. [C7] . [C12] . 
and Pi3,c given by Eq. [Aggf . and Eq. [61] from combining Eqs. [CTO] . [C12] . and Pi3,c. 

Figure IC] shows the dimensionless 3PT and linear power spectra, A2(fc) = k^ P{k)/ {2tt'^), for a matter component 
with pressure at different redshifts (z = 0.1, 1.0, 3.0, 5.0, 10 and 30) with kj = 1.0 and 3.0 h Mpc"^ The 3PT and 
linear power spectra are similar at the highest redshift, whereas the 3PT has significantly more power than the linear 
spectrum at larger wavenumbers as we go to lower redshifts. As a result, the filtering scale for a given linear filtering 
scale migrates toward larger wavenumbers in lower redshifts. 
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